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ABSTRACT 


An experimental investigation of typical foreign (impurity) defects 
and native defect generation by low energy protons has been performed* 
In- situ (vacuum) optical and spin resonance spectra were obtained 
which upon resolution revealed individual defect species and their 
densities. Analysis of these data provided a basis for correlating the 
growth kinetics of various color centers, including those associated 
with substitutional impurities and those associated with native lattice 
defects. In contrast to earlier theoretical work wherein discoloration 
was thought restricted to the proton range, it is now established that 
the dominant coloration mechanism consists of electron-hole pair dif- 
fusion into the crystal bulk (from a surface layer in which they are 
generated - proton range) followed by recombination and trapping in 
the bulk. This diffusion-dependence in the bulk (rather than proton 
displacement in the surface layer) greatly increases the influence of 
initial defect structure on the resultant optical damage. 
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1 . 0 INTRODUCTION 


The work reported herein represents an experimental investiga- 
tion of the discoloration of MgO by low energy proton irradiation* It 
constitutes a continuation of the first year l s theoretical study. Experi- 
mentation was restricted to MgO, both single crystals and polycrystalline 
material, as typical of a dielectric with application to thermal control 
coatings. 

Certain material preparation efforts were carried out. These 
included carefully controlled high temperature processing and selective 
doping of single crystal specimens with iron. 

Discoloration effects were introduced by bombarding specimen 
targets with 1. 0 keV protons (neutral beam) delivered at controlled 
rates in vacuum at room temperature. In- situ (vacuum) measurement 
of optical and EPR spectra were obtained and analyzed to demonstrate 
proton-induced defect species and generation rates. These kinetic 
data were then correlated in terms of individual color centers important 
to solar energy absorption. The relevance of accelerated laboratory 
simulation to real time in space was also indicated. 

Practical demonstration of the control of certain material 
parameters to obtain radiative stability is also presented. 


1 



2.0 EXPERIMENTAL 


The experimental portion of this effort was restricted to an 
investigation of damage induced in magnesium oxide (MgO) by low 
energy protons (1. 0 keV) incident at an order of 1000 x the solar wind 
proton flux. The effort involved measurement of changes in optical 
and EPR spectra resulting from controlled proton dose. 

2. 1 PROTON IRRADIATION FACILITY 

This facility is shown in Figure 1. Appearing in the lower center 
is an eight inch O.D. , horizontal, vacuum station connected to a 400 1/s 
ion pump (Hughes Aircraft Company, Vacuum Tube Products Division) 
at bottom center. A single Vac-Sorb, roughing pump is also shown 
attached to the station. 

Attached to the top of the station through a 1-1/2 O.D. entry is 
the proton source (seen in the center of the photograph). This is a 



Figure 1. 


Proton irradiation and optical 
measurement facility 
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trapped electron ionization source and was constructed by H. J. King 
and D.E. Schnelker of Hughes Research Laboratories. It is described 
in detail by King and Zuccaro' in their definitive study, "Solar Wind 
Simulation Techniques. " It provides a H + beam (>99 percent purity) 
of uniform intensity over a diameter of 2. 5 cm at the target specimen. 

Flux is variable over several orders of magnitude, and proton energy 
is variable up to 3. 0 keV. The associated gun electronics are shown 
at the right enclosed in a free standing console. These were designed 
and developed by J. Simpkins of Hughes Research Laboratories. They 
include a single feedback loop to provide flux control. 

At the upper left in Figure 1 is a Cary Model 14 spectrophotometer. 
Modifications of terminal optics in the instrument (enclosed in black 
painted housings) permit its scanning beams to be folded into an inte- 
grating sphere mounted within the vacuum station. 

Specimen movement within the station, e.g. between irradiation 
and optical measuring positions, is accomplished by means of a push- 
pull, rotary mechanism. Despite its use of a lubricated (Apiezon 
Number L) O- ring seal, this mechanism gave no evidence of introducing 
impurities which could serve to confuse the experimental work. 

Accessories included a Faraday cup (l.cm^ opening) movable to 
the sample position by a bellows- sealed transfer arm. A nitrogen gas 
aspirator was used to rough- pump the station, followed in turn by a 
sorption pump. Sample recovery in vacuum (at the end of an experi- 
ment) was accomplished by a 10 mm O. D. Suprasil II (fused quartz) 
tube attached to the station just beneath the sample position by means 
of a graded glass seal-Kovar- soft copper-2-3/4" O. D. , S. S. flanged com- 
pression port connection. Sample removal was effected by "pinch-off" 
at the copper tabulation. Manual control of the hydrogen gas flow to 
the ionizing gun was accomplished with a variable leak valve (Granville- 
Phillips Company Number 203019-02). 
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2.2 OPTICAL MEASUREMENTS 


A Cary recording spectrophotometer (Model 14) was used to 
obtain in- situ reflectance measurements before, during and upon 
termination of the particulate radiation. These measurements required 
about 3 0 minutes and were made on the specimen by quickly moving it 
from the irradiation position into the center of an integrating sphere 
located within the vacuum station. Accounting for the time required 
for instrument referencing, the specimen remained out of the p + beam 
for a total of about 15 minutes of the 30 minute measurement period. 
Because of noise associated with the e-source, it was found necessary 
to shut off the neutralizing e-beam during the entire 30 minute period. 

Certain terminal optics of this dual beam spectrophotometer were 
modified (see Figure 2) in order to permit its use for making reflect- 
ance measurements in-situ. First, the IR detector compartment was 
replaced with an IR lamp (6 v. , 18 a. ), thus reversing the IR beam 
path. Second, the UV/visible detector compartment was removed in 
order to permit both the reference and sample beam paths (ail modes: 
near UV , visible, and IR) to be folded into the integrating sphere 
(smoked MgO). This was accomplished as indicated by a pair of 
2 inch O.D. x 3. 84 F. L. spherical mirrors, a pair of 2 inch x 2 inch 
x 1/4 inch plane mirrors, and a single 2 inch O.D. x 5.91 F.L. 
spherical mirror. These mirrors were all first surface aluminum 
vapor deposited on optical grade pyrex (annealed), overcoated to pro- 
vide >85 percent reflectance between 0. 2 and 2. 5 p, and mounted in a 
light tight duct. The total increase in optical path was an order of 
30 inches. This increase, coupled with the necessity of beam stopping 
imposed by the dual requirements of reference/ sample beam segrega- 
tion and minimization of the sphere entry port, resulted in some signal 
attenuation. As a result, the lower limit of UV operation increased 
from a nominal 200 nm (without dry purging) to 235 nm. Operation 
in the visible region was unaffected. The upper limit of IR operation 
dropped to 1600 nm from a typical 2600 nm (2. 6 p) despite limited 
efforts involving the use of other IR sources of increased emission and 
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Figure 2. Schematic of modifications 


in Cary 14 optical system 


varying geometry. Since optically detectable defect structure in MgO 
in the 0. 8-2. 5p region is quite small, the spectral limitation in IR 
range did not prove a handicap. 

Interchangeable detectors were inserted into the detector port 
of the vacuum station. The pre-focused PbS cell used by the Cary 
for IR detection was retained. An EMI Number 9601B photomultiplier 
tube (S- 1 1 response) was used for near-UV/ visible detection, replacing 
the 1P28 (S-5 response) photomultiplier previously used. This change 
was designed to improve sensitivity in the visible region without 
serious loss in the near-UV region. 

Experimental single crystal specimens were arranged in an 
aluminum mount (0. 750 I. D. ) in a close fitting mosaic pattern of 
5-0. 5 cm x 0. 5 cm x 0. 09 cm (approximately) crystal wafers, resting 
on a second sur face, aluminized Suprasil mirror substrate. Upon 
termination, these crystals were quickly dropped into an EPR speci- 
men tube designed to sandwich all five together in a common <100> 
alignment. 

2.3 ELECTRON PARAMAGNETIC RESONANCE (EPR) 
MEASUREMENTS 

In the context of an experimental matrix of (a) selective impurity 
doping (b) controlled pre- irradiation processing and (c) control of 
absorbed radiation, EPR measurements can provide valuable data to 
supplement that obtained from optical reflectance spectroscopy. Such 
data prove helpful to mechanism studies in the following ways: 

a. Detection of paramagnetic centers 

b. Assist in identification of paramagnetic center 

c. Determine concentration of paramagnetic species and 
thus permit study in terms of generation or annealing 
kinetics 

d. Help to define the nature of the species surrounding the 
paramagnetic center. 
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Low temperature (e.g. 77°K) EPR measurements performed on 
irradiated specimens "in situ" represent the preferred measurement 
technique. Unfortunately the short range of 1 kev protons (the particu- 
late radiation source) requires a relatively large specimen surface 

area for exposure to radiation in order to induce a detectable number 

14 15 

of spins (e.g. , 10 to 10 spins). In turn, this imposes a need for 
a large-access microwave cavity and a special quartz dewar, both of 
which were acquired for this program. 

A modified Varian EPR spectrometer (model V-4500, reflection 
homodyne type) employing 100 kHz modulation of the magnetic field 
supplied by a 12 inch magnet (Varian V4012, variable yoke) was used. 
This system, including associated electronic gear, is shown in 
Figure 3. A block diagram of the spectrometer electronics used is 
given by Figure 4. Indicated is a swept frequency oscillator (SFO) 
used to determine the resonant microwave frequency for a given cavity 
loading at a given temperature (i. e. , specimen, specimen tube, dewar, 



Figure 3. EPR spectrometer system 
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Figure 4. Block diagram of EPR electronics 


cold gas, etc. ). Microwave power at this frequency is then introduced 
into the Varian microwave bridge, using as a source an ultra- stable 
microwave oscillator (Laboratory for Electronics Number 814-A-X-2 1 , 
8. 5 to lO.OkMHz) instead of the Varian klystron oscillator. The 
spectrometer is operated in the X-band near 9200 kMHz. Note that it 
was found necessary to alter the resonant frequency of the cavity (by 
the insertion of two aluminum rings) in order to locate the desirable 
TE Q13 mode near the mid-range of the LFE oscillator. A lock-in 
amplifier (Princeton Applied Research Number HR- 8) is employed 
for modulation and phase detection of the EPR signal which is dis- 
played on a strip recorder. A nuclear magnetic resonance precision 
gaussmeter (Harvey- Wells Number G-501)was used to measure the 
magnetic field intensity associated with the display of a given EPR 
signal. The proton probe of the gaussmeter is seen mounted (Figure 5) 
adjacent to the microwave cavity in the homogeneous field of the 12- inch 
magnet. This magnetic field value, coupled with the value of the 



Figure 5. Loaded microwave cavity 
arrangement in EPR spectrometer 
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microwave frequency, permits the determination of the g- values 
(see Appendix B) of the resonance lines in the EPR spectrum. 

Supported in the center of the magnetic field is a large sample 
access cavity (Varian V-453 5, . 2. 5 cm diameter access ) in which was 
inserted a quartz dewar (Varian P/N 909005-01, 11. mml.D. ). Cold 
nitrogen gas (~80°K), from a liquid N 2 chiller- transfer system, was 
passed up through the dewar and around the evacuated specimen tube 
contained therein. This arrangement of specimen tube-dewar- 
microwave cavity - cold gas system is shown between magnet pole 
pieces in Figure 5. The Q of this cavity- dewar- sample tube configu- 
ration has been measured to be about 6600 and provides an optimum 

14 

sensitivity of approximately 5x10 spins. 

2.4 MATERIALS: EFFECTS AND PREPARATION 

As a preface to this section it is important to point out the 
dominant role played by lattice- soluble impurities in determining the 
optical changes introduced by radiation. Further, as will be demon- 
strated in Section 3. 0, the magnitudes of these effects may be signifi- 
cantly influenced by the thermal treatment imposed during or subse- 
quent to initial preparation. It is equally important to recognize that, 
although MgO has received more research attention than any other 
alkaline earth- or refractory- oxide in terms of radiation- induced 
optical changes, precise definition of the effect of individual impuri- 
ties (either substitutional or interstitial) in varying matrices of other forign 

(extrinsic) and native (instrinsic) defects remains a distant goal. Typical of 

2 

the problem, as recent as 1967 Chen and Sibley reported that 
center generation in MgO by ionizing radiation (important here) was 
crystal- dependent. The nature of this dependency could not be identi- 
fied. They did find that very fast quenching (e.g. , 23°C/sec) from 
temperatures between 500 and 1250°C prior to irradiation could 
increase resultant V^ densities by as much as five-fold, depending on 
the source of the crystal. They also reported that Co, Cr, Mn, and 

3 

Ti act to suppress V^ center production. Dr. J. L. Kolopus , an 
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associate in the above work, states this particular problem has not yet 
been resolved and probably will require an uLtra-pure' (e.g. , 1-2 ppma 
total impurities) target material for its resolution* Parallel work 
reported by Glass and Searle**, while not conclusive, did indicate that 
an increase in isolated magnesium ion vacancies results from such 
thermal treatment, yielding an increased density of centers. They 
pointed out that u cation vacancies and impurities are essentially 
immobile in MgO below about 400°C* 11 

Despite the above mentioned data limitations, it is believed that 
the present state-of-the-art in MgO is sufficiently advanced to permit 
significant progress toward obtaining a truly radiation- stable pigment. 
The relatively simple, cubic symmetry of MgO provides another 
reason for selecting this compound, from among other white dielectric 
oxides, for initial experimentation* The work of the program was thus 
confined to MgO, powder and single crystal. 

Reagent grade powder {Mallinckrodt Chemical) of relatively low 
purity was hydraulically compressed (50-70 kpsi) at room temperature* 
Easy delamination of the pressed discs permitted random chip selec- 
tion from within the pressed body. A weighed amount of these irregu- 
lar chips ('-O. 4 cm x 0. 4 cm x 0. 13 cm) constituted the proton target. 
They were not annealed or given other thermal treatment prior to 
irradiation. The intent was to use the impure powder spectra (both 
optical and EPR) for initial study as a guide to the overall defect 
structure. Despite the availability of high purity powder, subsequent 
experimentation was carried out with single crystal MgO in order to 
obtain n cleaner M spectra. This is especially true of EPR spectra 
because of orientation dependency associated with anisotropic defects. 

Single crystal MgO was obtained from two sources. High and 
moderate purity material was provided through the courtesy of C.T. 
Butler from the Research Materials Program of the Solid State Division 
of the Oak Ridge National Laboratory. Low purity material was 
obtained from the Norton Company (optical quality 11 Magnorite 11 ). In 
terms of possible substitutional replacement -of lattice ions, the impuri- 
ties selected in Table 1 are likely candidates. In this manner they can 
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Table 1. Concentration of Selected Impurities in Various MgO Specimens 

(Parts Per Million Atomic ) 


Element 

Detection 

Limit, 

ppma 

Compressed 
Powder- Very 
Low Purity 
(-43) 

ppma (UP cit?) 

Single Crystal, ppma 

Low Purity- 
N-? Processed 
" (-66) 

Moderate Purity: 
Air-Quenched 
f71) 

High Purity- 
As Received 
(-55) 

High Purity- 
No Processed 
“ (-63) . 

High Purity: 
Fe- doped, 
Air-Quenched 
(-67) 

El 

1. 


4. 9 

5. 1 

6 

6. 

6- 

B 

19. 

nd 

nd 

nd 

nd 

nd 

18. 6 





* 

* 

* 


N 

-- 


## 

2.9 

7.2 

7.2 

7.2 

F 

-- 


O* * 

<1. " 

<1. * 

<1. * 

<1. V 

Na 

0. 1 , 

83 0. 

28. 1 ' 

2.2 

9.3 

12. 2 

nd 

Al 

3.7 

12.7 

16.4 

11.2 

n. 9 

nd 

10.4 

Si . 

14. 4 

108. 

20. 8 

-- 

23. 6 

38.7 

36. 5 

P 

-- 


## 

2. 6* 

26. " 

26. " 

❖ 

26. 

Ti 

4.2 

nd 

1 3* 

1.3" 

nd 

nd 

nd 

V 

4.0 

nd 

nd 

nd 

nd 

nd 

nd 

Cr 

0. 4 

nd 

1. 3 

nd 

nd 

nd 

nd 

Mn 

1. 5 

i 

nd 

20. 6 

nd 

nd 

21. 6 

nd 

Fe 

1. 8 

7. 6 

39. 7 

8 0 

2. 0 

5.4 

541. 

‘ Co 

1.7 

nd 

nd 

nd 

nd 

nd 

nd 

Ni 

1. 4 

nd 

1. 6 

nd 

nd 

nd 

nd 

Cu 

0. 2 

1.2 

0. 3 

0. 4 

0.3 

0. 4 

0.4 


nd - not detected (below detection limit) 
^Estimate 
^Not determined 



generate intrinsic defect structure initially and during irradiation 
quench- in radiation- induced defect structure as trapping sites, indi- 
vidually or in association with intrinsic defects. With the exception of 
Li and Na {determined by flame photometry and atomic absorption 
analysis respectively) and N and F (estimated from typical analyses on 
this material previously made), the remaining elements were obtained 
from arc emission analyses* It was originally intended to use spark 
source mass spectrometry in order to obtain <1 ppm detection limits, 
but as the work proceeded such improved sensitivity did not seem 
required, with the possible exception of boron* 

2,4*1 Hydrogen 

Hydrogen analyses were not obtained. This may represent an 
important omission from Table 1. The degree of importance probably 
relates to the density of OH , either resulting from (a) original prepar- 
ation procedures or from (b) proton irradiation. The latter case might 
include both proton (H + ) stopping in the "surface 11 : 


H + + o 2 "— OH 


and diffusion of charge equilibrated hydrogen atom (H -e) interstitials 

into the bulk with capture at a V, center, i. e. , the latter described as 
• . n 


a hole h localized on an O adjacent to a magnesium ion vacancy, V 


Mg' 


° 2_ h- V, , 
Mg 


, T ,+ 1 

+ H -e 


OH -V 


Mg 


4 

Henderson and Wertz review this area in detail because "it has 
led to a description of the geometry of one stable site for hydrogen in 
an oxide. 11 Current infrared studies indicate that substitutional OH in 
isolated sites is displayed optically by a broad band centered at a wave 
number of 3400 cm ^ (\«2. 94(r), with the hydrogen atom suspected of 


14 



being directed inter stitially . Additionally, in unirradiated crystals 
containing OH a strong, narrow band is displayed at 3296 cm - * 

(X«3. 03p). This has been associated with the axially aligned "trio", 

2- + " ii 

O H (equivalent to OH )• Upon ionizing irradiation, the 

3296 cm~^ band is shifted to 3323 cm“l (3. Op) with attendant generation 
of a strong EPR signal. This latter is associated with the axially 
aligned "quartet", OH -V^^'-O , where O” represents a h* localized on 
a normal lattice oxygen ion. Note that all three above spectra are 
more prominent in MgO single crystal prepared from Mg(OH), powder 
than from MgC0 3 powder. The aniosotropic defect, OH'-Vj^-O , 
has been termed the Vq.^ center. 

Interestingly enough, examination of Chen and Sibley's^ data 
(Figure 6) reveals that ionization - induced center saturation is 
inversely proportional to OH density prior to irradiation. Their data, 
which include crystals from various sources with random variations in 
total iron of 5-100 ppm and total chromium of 1-20 ppm etc. , are 
presented in Figure 6. If the ultimate fate of a significant portion of 
low energy protons stopped in a target is to form OH , this may also 
explain part of the difficulty during this work of generating larger 
densities of Vj centers by proton irradiation. 

The role of OH , both isolated and associated with V" , in 

Mg 

radiative coloration of MgO warrants further attention. Such investi- 
gation was not pursued during this work because other aspects 
appeared more pressing. 


2.4.2 Elements of First Transition Series ' 

Ions of this series have received the most attention and have been 
shown to significantly affect radiation coloration processes. The more 


important of these are as follows: 
Ion 

c 

Preferred Valence State 

Fe 1+ , Fe 2+ , Fe 3 + 

Fe 2 + 

Cr 2 +, Cr 3 + 

Cr 34 

Ti 3+ , Ti 4 + 

Ti4 + 

Mn 2 + 

Mn 24 

Co 1+ , Co 2 + 

Co 2 + 
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MAX. ABSORPTION COEFFICIENT FOR Vj^ CENTER, a MAX IN CM' 


r-H 



0 0.5 1.0 1.5 2.0 2.5 . 3.0 


ESTIMATED MAX. ABSORPTION COEFFIC IENT FOR OH' BAND, a QH IN CM " 1 

7 

Figure 6. Chen and Sibley data for y -radiation induced maximum Vi center absorption 
as a function OH" content in MgO single crystals from various sources 



There is limited substitionaL solubility of V 3+ /V 2+ (not of V 5+ ) 6 but 

no indication of being of real importance to radiation induced defect 

kinetics. The divalent ions o£ Ni and Cu also appear to have limited 

solubility in the MgO lattice but no apparent importance to radiative 

coloration. Note that the preferred valence state of the iron and 

titanium are subject to thermal history. Thus for the region of impur- 

ity levels shown in Table 1, oxidation at 1200-1500°C favors the states 
3+ . 4 "P 

Fe and Ti , while reduction (in vacuum, inert gas, or hydrogen) at 
o Z + 3 + 

1200-1500 C favors Fe and Ti . In terms of iron , such partition of 
valence states by means of thermal treatment and environment control 
is illustrated by the following experimental data (Table 2) on both low 
and high purity single crystal MgO. 
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Table 2. Processing Effects on Fe /Fe Distribution in 

Single Crystal MgO 


Process 

Ions * 

7 3 

/cm 

v 3 + 

F e 

_ 2 + 
F e 

Low Purity Crystal 



As received (expt. -50) 

4. 8 x 10 18 

0. 0 

Soak at - 145 0°C in N 2 for 4 hrs. 
followed by intermediate cooling 
rate ^ 2 5°C / min (expt. -66) 

1 O 

0. 02 X 10 

, - ir ,18 
4.2 x 10 

High Purity Crystal 



As received (expt. -55) 

2. 0 x 10 17 

0.7 x 10 17 

Soak at ~1450°C in N 2 for 4 hrs. 
followed by intermediate cooling 
rate = 25°C/min. (expt. -63) 

17 

0.4 x 10*' 

17 

5.4 x 10 x ' 

* 14- 

Fe not considered present. The variation in total iron content 

(Fe^ 1 * + Fe 31 *) observed for a given materiaL results from (a) con- 
lamination during processing, (b) errors in instrumental analysis, 
and (c) concentration variations within a single ingot. 
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The ions Cr 
such crystal treatment 


Mn^, and are typically not affected by 


2,4.3 Alkali Metal Cations 

£ 

Adam and Brindley investigated a range of monovalent, divalent, 

trivalent, tetravalent, and penta valent cation and cation pair solutions 

1 -J* 3 4 

in crystalline MgO. The Li -A1 pair was reasonably soluble but 

, t 1 • + -j-. 14-.34 14 34 T*^ r ^ j 

the Na -A1 pair was not* Li -Cr , Li -Mn , Li -Fe , and 

1+ 4+ 1+ 54 1+ 44 

Li -Ti pairs were also quite soluble. The Li -V and Li -Mn 

pairs were not. From their data it seems reasonable to assume that, 

of the alkali metal ions, onLy Li provides cation charge compensation 

in the lattice for substitutional ions of valence greater than 24. The 

importance of such charge compensation may be understood by recog- 

- . 3 4 

nizing that for each two trivalent cations (e. g, , Fe ) which are not 

IT 

compensated a is generated. This is a potential hole trap. 


2.4.4 Miscellaneous Cations 

3+ 4 -j- 3+ 3 + 

, Lattice effects of A1 , Si , P , and B , either in substitu- 
tional or interstitial positions, are not yet sufficiently understood to 

gage their role in -r egulating radiative coloration processes. However, 

34 

it is not unreasonable to suspect that A1 , for example, might par- 
ticipate as an* electron trap or as a recombination center. 


2.4.5 Anion Impurities 


Trivalent, substitutional anions (e.g. , N ) can generate 

1- U 

while monovalent anions (e.g. , F ) can oppose such oxygen vacancy 
generation. In addition to generation of trapping sites, each may 

participate directly as a defect state. Thus the V center has been 

4 " * 

reported to consist of: F -O . . This center is observed to 

be the analog of the center. 

■Nothing is known of defect states involving nitrogen ions. 
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2.4. 6 Fe -Doping Procedure 

3 + 

Selective doping of Fe into a high purity MgO lattice offers 
the potential for: 

2 + 

a. An increase in the concentration of Mg vacancies 
provided it is not accompanied by an increase in Li 
concentration. 

b. Controlled introduction of a recombination center and 

a carrier trap at the same time. The nature of the 

trap appears to be a function of the level of the 

doping, net electron trapping being associated with 

[Fe^ + ] > - 4 x 10^/cm^ ( Net hole trapping 

° 3+ 

seems to occur below this Fe level with the impurity 
matrices shown in Table 1. 

A thermal diffusion procedure was developed for this doping. 

The first step consisted of thermal vaporization of high purity iron 

- 5 

(6-9’s purity) at 10 torr on to freshly cleaved specimen wafer s of 
high purity MgO (0. 75-1 mm thick). Coating thicknesses of the order 
of 2500 A were used for high level doping. Lower iron doping levels 
would require commens urate ly thinner coatings. These crystals were 
stored in a dessicator until packaged for use just prior to thermal 
tr eatment. 

Packaging for thermal treatment was designed to minimize 
impurity contamination from the high temperature ambient. Working 
in a n clean cabinet 11 the specimen crystal wafer s were first enclosed 
in a platinum foil capsule measuring 3/8 inch O.D. x 1 inch 
(0. 0005 inch foil thickness). This first capsule was then enclosed in 
a second slightly larger Pt foil capsule. This package was next 
immersed in high purity MgO powder contained in a 3/4 inch OD. x 
1-5/8 inch Pt container and then covered. This final package was 
loaded into a MgO crucible (Leco Number 528-7 01, 99. 5 percent wt. 
MgO typical, 1-1/4 inch O.D. x 1-3/4 inch) fired to 1500°C prior to 
use. 


; [ W' 
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Immediately before carrying out the diffusion experiment, the 
loaded MgO crucible was transferred from a desiccator to a graphite 
susceptor (baked out at 110°C for 24 hours), covered with a graphite 
lid containing a pyrometer view port, and placed in an induction furnace 
with rotary pedestal (Future craft). The system was closed and care- 
fully purged with dried N 2 gas. Continuing this inert gas purge, the 
temperature of the Pt package was raised to ~T460°C in approximately 
one hour. A temperature of 1460-1480°C was maintained for nine 
hours, at which time the crystal package was rapidly cooled to <500°C 
at an initial rate of ~6°C/ sec which then decreased to ~2. 5°C/sec. 

The cooled, loaded MgO crucible was quickly transferred from the 
induction furnace to a desiccator and taken to a clean cabinet. Here, 
the entire MgO crucible was wrapped in 0. 5 mil Pt foil. 

This latter package was next placed in a globar furnace (SiC ele- 
ments) for selective oxidation of the iron introduced into the lattice. 
Employing an air ambient, the temperature was raised to ' vj 1480°C 
wi.thin about five hours. After soaking at this temperature for four 
hours, the hot MgO crucible was quickly removed from the furnace 
and air-quenched on a heavy aluminum plate (cooling rate ^10°C/sec 
initially, dropping to ~-2. 5°C/sec at 500°C, 

Table 3 indicates the effectiveness of this doping procedure 
and is based on optical, EPR and spectrographic analyses. 
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Table 3. Fe -Doping of High Purity Single Crystal MgO 


Process 

* 3 

Ions /cm 

TT 3 + 

F e 

TT 24 

F e 1 

Undoped (-55) 
Fe^ + -Doped (-67) 

2. x 10 17 
1. 5 x 10 19 

17 

0.7x10 ' . 

4. 3 x 10 i9 

. — ■ 

^ 14 

Fe not considered 
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The Fe content is noted to have been increased two orders of 

magnitude* Reference to appropriate experiment numbers (in 

parentheses) in Table 1 indicates that this increase was accomplished 

3 + 

with little significant contamination. An increase in the [Fe ] rela- 

2 + 

tive to the [Fe ] in the doped specimen can be effected by a 
longer soaking period during oxidation. 

Uniformity of bulk distribution was not determined* This is an 
important consideration. However, careful visual examination of the 
doped crystals revealed a uniform yellow coloration, excellent trans- 
parency (except for a minor number of striae initially present), no 
cloudiness, and no obvious clustering or other irregularities. While 
not conclusive, such cursory examination seemed to preclude the need 
during this phase of the effort for more sophisticated examination of 
doping uniformity, e* g. , by means of direct imaging mass analysis. 
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3. 0 COLOR CENTER GENERATION IN MgO 


In addition to several preliminary experiments with MgO powder 

to permit check-out of irradiation and detection facilities, a total of 

seven, radiative coloration experiments were performed. Each was 

continuous, ranging in duration from 6 to 12 days. All were conducted 

*8 ~5 

at room temperature in a vacuum of ^10 torr backfilled to 10 
torr with high purity gas (Matheson Co. , 99. 999 5 percent vol. ). 
Proton energy was maintained at 1. 0 keV in all cases. Proton flux 
was varied over the range 3x10^ to 4 x 10^ p*/cm 2 -sec. 

Depending on which value is accepted for the solar wind proton flux, 
experimental proton irradiations were carried out at rates of the order 
of 1000-1500 X solar. In parallel, an~l. eV electron (e 1 ) beam, (non- 
luminous in the visible wavelengths) was directed at the target at 
an intensity slightly greater than that of the proton beam. Flux values 
were controlled within the limits of ~il0 percent. Reflectance changes 
were measured as often as once /day or as otherwise deemed necessary 
to monitor the color changes produced. Upon termination of an experi- 
ment, corresponding reflectance and EPR spectra were obtained. 

3. 1 POWDER EXPERIMENTS 

The first experiment was conducted with powder. An impure 
material (reagent grade, see Table 1) was selected deliberately and 
compressed into chips, without subsequent thermal annealing, to 
further increase the defect structure. Irradiation of such material 
was expected to provide a broad spectral look at the distribution of 
impurities one might encounter under worst conditions. Such 'exami- 
nation might then suggest means of keeping the problem tractable and 
yet confined only to those areas pertinent to solar thermal control. 

The absorption change generated by a total proton fluence of 

17 4. 2 

2. 6 x 10 p /cm is shown in Figure 7. Estimated Gaussian resolu- 
tion (see Appendix A) of the overall spectrum into component spectra 
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PHOTON ENERGY, eV 


200 250 300 350 400 500 600 700 800 

WAVELENGTH, nm 

Figure 7. Absorption change in MgO powder (compressed) induced by a 
1. 0 keV neutralized proton beam at 20°C and 10“ torr H£ (F1799-43) 




+ 34 * 

is indicated. The F -center and Cr * spectra are only approximate. 

The importance of this experiment was to point out the following: 

a. The major portion of the coloration even in a very impure 
dielectric, is displayed in the near-UV. 

b. Component resolution at photon energies (hv) greater than 

~5. 0 eV is seriously handicapped (if not impossible) at 

this time because of the lack of model information of 

resonant color centers in this region. Henderson, King, 

7 

and Stoneham have made tentative suggestions in this 
regard. 

c. A number of unidentified spectra are displayed at the photon 
energies corresponding to the visible wavelength region. 

These overlap the identified center. Since the 
resolution of the band is important to the attempt 

to study defect kinetics, reduction of adjacent overlapping 
band structure might be attempted by resort to higher purity 
materials. Additionally, since EPR powder spectra are 
more complex than those for the single crystal, attention 
was further directed toward higher purity single crystals 
for subsequent experiments. 

d. As far as near-infrared defect display is concerned, little 
was observed. Many workers have reported a band at 

1. 2 eV (1. 03fJi), and Henderson and Wertz have noted that it 
changes in concert with the 3. 6 eV band. This band was 
not observed in this work. Small initial bands at ~1. 55 eV 
(0. 8p) and ~0. 9 eV (1. 38p) were observed regularly, but were 
unchanged by irradiation. These are likely associated with 
free water and/or CO^- 

3. 2 SINGLE CRYSTAL EXPERIMENTS 

3. 2. 1 Effect of Controlled Chemical Processing 

A primary consideration to color center experimentation is the 
defect structure {both intrinsic and extrinsic) resulting from initial 
preparation and as modified by post-processing. This matter has 
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been discussed earlier but is such a dominant consideration as to 
warrant both restatement and illustration. Present limitations in 
defect modeling and precise characterization techniques limit the 
guidance available to preparatory efforts. Thus, (a) careful planning 
within the context of available data and (b) very careful handling 
procedures are an absolute requirement. Illustrative of item (a) is the 
previous discussion of the possible effect of OH concentration on the 
radiation-induced density of the center which absorbs solar energy 
near the peak of the solar spectrum. The presence of OH~ appears to 
suppress Vj center formation by favoring center formation, the 

latter absorbing optical energy is a narrow band centered at 3. Op' 

{i. e. outside the solar region). Recognizing that other factors might 
also enter, at least the option of introducing OH" (say by starting with 
purified MgfOH^) in order to suppress radiative' formation of Vj centers 
should be considered prior to material' preparation. 

Care in material handling, the second item above, as well as the 
option of post-preparatory processing are well illustrated' by Figures 8, 
9, and 10. The crystals used for the two experiments involved (-55, 
and -63) were' all cleaved from a single, small section within one high 
purity ingot. Careful reduction (and annealing), as indicated-partly by 
little increase in impurities in Table 1, is demonstrated by comparison 
of Figures 9 and 10. Note that the processed material was -subjected to 
over twice the proton fluence imposed on the as-received material. 
Despite this, radiative-induced optical absorption in the near-UV is 
reduced by a factor ranging from 1/3 to 1/10. Thus, chemical reduc- 
tion of the iron 3+ state, without attendant contamination or other 
deleterious effects, contributes by decreasing radiative -coloration over 
broad regions centered at 4. 25 eV and 5. 7 eV (the latter Fe band not 

being shown). The 3. 6 eV band, although not yet identified, is 

8 

considered by Hansler and Segelken to represent absorption by a center 
consisting of an electron trapped at an intrinsic lattice defect. Com- 
parison of Figures 9 and 10 would indicate that this defect was effectively 
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Figure 8. Change in optical absorption induced in processed , high purity MgO single 
crystal by neutralized 1. 0 keV proton beam at 20°C and 10“^ torr H £ (F1799-63) 
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Figure 9. Change in spectral absorption coefficient induced in high purity MgO single 
crystal by neutralized 1.0 keV proijon beam at Z0°C and 10"^ torr H 2 (F1799~63). 
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Figure 10. Change in spectral absorption coefficient in,duced in high purity MgO single 


crystal -by neutralized 1. 0 keV proton beam'at 20°C and 10" 5 torr {Fl799:-55) 



removed by annealing as a result of processing at 1450-1500°C in 

nitrogen gas. Other, unidentified near-UV bands, either generated or 

which grow during p -irridation, are seen to be similarly reduced as a 

result of such controlled chemical processing. 

Similar decrease in spectral absorption in the visible region is 

also noted. The center which requires a density of Vj^ was not 

observed in the processed crystals. This fact indicates that the 

[ V 1 * , i.e, prior to irradiation, was negligible and results from a 

M gJo f 3+1 

minimized Fe I and unidentified annealing effects. The effect of 
p^-generated V M ^g in the proton range is discussed later. 

As mentioned earlier, negligible radiative coloration was 
observed in the near-IR. This was true even in the case of very 
impure, highly stressed materials. 

Please note that the absorption coefficients which are referred 
to in the above figures are not true absorption coefficients. That is, 
they do not represent a uniform coloration of the bulk thickness. This 
results from proton stopping in the surface layer and subsequent 
charge carrier diffusion into the bulk to some undetermined distance 
less than the specimen thickness. Thus the A a ^values represent a 
mean value over the total crystal thickness. Since crystal thicknesses 
varied among experiments only ±10%, such values provide a useful 
basis for comparing defect spectra. 

3.2.2 Effect of Selected Impurities 

3.2. 2. 1 Iron. Comparison of the crystals used in experiments -63 and 
-66 on an impurity basis only (see Table 1) reveals one major variation, 
namely the total iron content. Both were processed simultaneously ; 
after which the [total Fe] ~ p [total Fe] Processing was designed 

to minimize the initial density of the 3+ state of iron and anneal out 
some of the intrinsic defect structure. Comparison of the optical 
absorption spectra for the unirradiated specimens (see Figures 8 
and 11) indicates considerable similarity. Only after comparable 

17 } 2 

proton irradianc'e (>10 p /cm ) does the difference between the 
relatively pure (low iron content) and the impure (high iron content) 
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IRRADIATED 

(1.4X10 17 p + /CM 2 ) 


UNIRRADIATED 


SINGLE CRYSTAL: NORTON CO. 
"MAGNOR1TE 11 PROCESSED IN N ? AT 
1450-1500 °C; x -0.096 CM 
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WAVELENGTH, nm -► 

Figure 11. Change in optical absorption induced in low purity MgO single crystal 
by a neutralized 1. 0 keV proton beam at 20^0 and 10~5 torr H 2 (FI 7 99 -66) 







crystals manifest itself. Thus proton irradiation increases the solar 

9 

absorptance Og (see Blair 7 for spectral weighting factors) of the pure 

crystal an order of <5% from ~0, 03 to ~0. 031 . Compare this to an 

increase of the order 200% for the impure crystal, from ~0.025 to 

~0.075, for an even lower proton fluence. 

The reasons for this divergence may be observed by comparing 

component defect spectra (see Appendix: A for basis of Gaussian 

resolution technique) shown in Figures 9 and 12, In the low purity 

crystal, certain features are obvious as follows: 

34- 

a. Growth of the Fe bands (i.e. resulting from net h'-trapping 

2+ 

by Fe recombination centers), 

b. Filling of intrinsic, electron traps, e, g, by and yielding 
absorption near 5, 0 eV (F and F centers) and by absorption 
centered at 3, 6 eV, 

c. Filling of intrinsic, hole traps, e. g, by yielding 

absorption centered at 2,3 eV (V^ center). 

In the xase of the high purity crystal, all of the above are suppressed 
to varying degree. 


Selective doping with iron to a relatively high density can also 
suppress radiation-induced defect density. Regardless of processing, 
both and |Fe^ + J Q are so large that, acting as recombination 

centers, they are able to suppress trapping of proton-generated charge 
carriers, both e J and h\ Thus proton-induced coloration is small. 

Such high level Fe -doping is indicated in Table 1, noting that the total 
iron content of crystals in experiment (-67) was increased to~100 X 
that of the crystals in experiment (-55), other impurity concentrations 
remaining about the same. The penalty associated with such bulk- 
suppression, in terms of thermal control, is very broad, initial,, UV 

Or 

absorption tailing strongly into the visible region. This Fe - related 
absorption would not prove acceptable in a thermal control application. 

19 ' o 

Because at densities of 5.8 x 10 total Fe/ cm (= 541 ppma), the 
34- 2+ 

Fe /Fe densities dominate the electronic processes, this material 
did not prove too useful in studying defect kinetics. In concluding, note 
that Haxby^ also observed this effect of iron-doping (0, 01 to 0. 40% wt. ) 
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CHANGE IN ABSORPTION COEFFICIENT Aa , AT 20PC, CM 






in regard to color center suppression (in the bulk) during X -irradiation 
of MgO single crystals (compare his Figures 3 and 10). 

3. 2. 2.2 Hydrogen as OH”. A broad-, near-IR band centered near a 
•wave number of 3400 cm" ( \= 2. 94p) has been associated^ with 
substitutional OH in isolated sites. Glass and Searle^ have also 
attributed a nearby band (at 3700 cm ) to Mg(OH)^ precipitates in MgO. 

Optically resonant associated (or aggregate) centers formed by 
adjacent - OH pairs with an adjacent trapped h* have been 

discussed earlier in Section 2. 4. 1. They are discussed further in 
Section 3. 3 in terms of the EPR work. Neither these aggregate centers 
nor the isolated OH center were examined by optical spectroscopy 
during this work. However, they appear to be of importance. 

Two important conclusions derive from the above discussed 
data as follows: 

a. The term purity must be considered in a very selective 
manner and is only relevant to the problem of radiative 
coloration when employed in the context of a given crystal 
lattice. 

b. Despite current gaps in color center knowledge, this body 
of information can be introduced into preparatory .and 
(post-processing procedures in order to suppress radiative 
coloration of a dielectric, e.g. MgO. Preliminary applica- 
tion of such understanding toward development of a radiation- 
stable pigment for spacecraft thermal control is illustrated 
quite well by Figure 8. 

.3. 3 POINT DEFECT IDENTIFICATION BY SPIN RESONANCE 

Electron paramagnetic resonance (EPR) measurements were 
carried out on the crystal specimens to provide supplementary defect 
identification and measurement for the optical data. EPR measure- 
ments were only possible after experiment termination. They were 
arranged so as to prevent post-experiment optical and/or thermal 
bleaching. Thus EPR measurements and terminal optical. spectra are 
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considered to correspond^ both being made in- situ (vacuum) . Handicapping 
these EPR measurements were such factors as the following: 

+ 

a. Microwave power saturation of selective defects, e.g. F . 

b. Orientation problems associated with various anisotropic 

3 + 

centers, e.g. Fe and the various V -centers. 

c. Errors introduced into defect density calculations because 
of poorly defined spectra, especially spectral tails. 

d. The presence of defects with overlapping spectra, making 
resolution of components difficult and subject to error. 

e. Inability to display certain potentially observable valence 
states satisfactorily except at 4°K (e.g. Fe^) and others 
not at all (e* g. Fe^*). 

These and other problems serve to limit the potential usefulness of the 
EPR spectra. Despite such shortcomings, these data did serve as a 
valuable supplement to the optical measurements. 

Defect densities calculated from typical EPR measurements are 
compared below with results derived from corresponding optical 
measurements (see Table 4). 


Table 4. Selected Defect Determinations Obtained From 
Corresponding EPR and Optical Data 



17 3 

Defect Density, 10 /cm 

Fe 

_ 2+ u 
Fe * 

V 1 

Cr 3+ 

Experiment 

(-50) 




t 

Optical j 

Unirrad. 

48.0 

0.0 

0. 0 (est.) 

nd 

p^-irrad. 

34.0 

14. 0 

0.2 

nd 

EPR 

Unirrad, 

54. 0 

0 . 0 

nd 

13. 0 

p + -irrad. 

31. 0 

23.0 

nd 

11.0 

Experiment 

(-71) 





Optical 

[ Unirrad. 

4. 6 

4. 1 

0.0 (est.) 

nd 

| p^-irrad. 

is. 4 : 

2.3 

0.2 

nd 


Unirrad. 

_ « 



nd 

3. 

EPR 

p + -irrad. 

! 

-- 

nd 

3. 


y-irrad. (satn.) 

0.7 

8. 0 

0.8 

0.7 


-Obtained by subtracting Fe 3+ from total iron density, 
nd - not detected 


35 




As indicated by these typical data, EPR and optical results were 
generally consistent. 

Gamma radiation (in vacuum) of corresponding crystal specimens 
were frequently carried out to examine bulk saturation effects. Typical 
dose was 15 Mrad of 1 MeV Y-radiation. The Y-induced EPR spectrum 
may be compared in Figure 13 with the spectrum of the unirradiated 
specimen used in experiment (-71). The region covered ranges from a 
low field side of H = 3254. 94 gauss (g = 2. 0383) to a high field side of 
3351.40 gauss (g = 1.9796) in the direction shown. The following 
defect density values attained at saturation by ionizing radiation were 
estimated from component V^, V^, and spectra in Figure 13 as 
follows : 


~ 3. 53 x 10^ /cm^ 
~ 0.79x 10^/cm^ 
- 0. 32 x 10^/cm^ 


Total » 4. 6 x 10^ /cm^ 

These data provide information pertinent to defect studies in the 
following manner: 

a. Comparison between generated by p + - and by 

Y-ir radiation indicates the approach to saturation in 
regard to this center, i.e. 25 percent in this case. 

b. Based only on impurity -introduced cation vacancies, 
one may estimate jV,, -1 as an order of 10 x 10 ^/cm^, 

2+ - L M SJo 

Mg vacancies introduced thermally (as during 
preparation) would add to this value of course. In any 
event, one would expect a significant portion of these 


[ v oh] 

M 

O] 
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MEASURED SPECTRUM 

RESOLVED COMPONENTS 


Figure 13. Effect of gamma -irradiation (in vacuum at 20°C) 
on EPR spectra (79°K) of air-quenched, 
moderate purity, single crystal MgO 




vacancies to become converted to V-centers as a 

consequence of K-trapping during Y-irradiation. In 

this regard compare the saturation density of V-centers 

17 3 

(measured by EPR above - 4. 6 x 10 /cm )' with the 

2 + 173 * 

Mg vacancy density of >M0. 0x10 /cm . 

By means of the Smakula -Dexter equation relating defect 

density with optical resonance, one may estimate 


V 'max* Separate estimation of allows comparison 


with the data in Figure 6, Thus: 


N - S 


a 


max 


where 


N = defect density, no, /cm 0 

1 6 > 2 

S = Smakula -D exte r band constant = 6, 4 x 10 cm 


for center 

a max = absorption constant at band peak, cm 


-1 


From above, using data presented earlier and in Appendix A: 


(\) 


0. 79 X 10 17 _ . -1 

= 1.2 cm 


max 6.4 x 10 


16 


Approximate ^OH“J in MgO by an expression given by Kats^ 7 for £oH"^ji 
quartz crystal as follows: 


N oh - = 6. 32 x 10 


17 


(“oh) ; 


max 
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[ 1 17 3 

Vqh centers ~5 x 3. 53 x 10 /cm # 


["oh] 


5 x 3. 53 x 10 


17 


max , ln 17 

6. 32 x 10 


~ 2.7 cm 


-1 


These corresponding maximum values of and Gqjj ar e plotted in 

Figure 6 and show reasonable agreement. 

In EPR work, the detection threshold is typically a direct 

function of absorption signal half width. The detection threshold for 

14 

the EPR system used is estimated as 5 x 10 spins. On the basis 

3 

of the typical sample volume « 0. 1 cm , a detection threshold of an 
15 3 

order of 5 x 10 defects /cm , MgO is thus established for defects 

3+ 

displaying relatively narrow lines, e.g, Cr in cubic symmetry, 

____ ^ 

It was found that broader lines, e. g, Fe , could not be observed 

1 6 34 3 

below an order of 5 x 10 Fe /cm MgO. 
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4. 0 COLOR CENTER GENERATION KINETICS IN MgO 


1 2 

Theoretical studies performed during the first year of this 
effort indicated that defect kinetics associated with low energy protons 
(in dielectric targets of thickness » proton range Rp) were regulated 
by proton- generated lattice displacement within a "surface" layer 
defined by Rp. Thus in the case of MgO, the center density > 

deemed important because of its optical resonance in the visible spec- 
trum, was indicated to track the density of magnesium ion vacancies 

primarily determined by time- integrated proton dose. Steady 
state or saturation defect densities were expected when high probability, 
V Mg ” Mg'i pair separation was .attained as a result of irradiation. 

Estimates of such Frenkel pair separation led to expected, saturation, 

21 3 

defect densities approaching 2. 7 x 10 /cm as a limit. It was appre- 
ciated that charge exchange processes involving impurities were impor- 
tant in terms of their function as individual color centers. In the case 
of dielectrics these centers typically display near-UV absorption and 
are thus of minor importance to sola'r absorption. It was believed the 
primary role of substitionally soluble impurities was to participate in 
regulating the available density of holes (h) or electrons (e -)-at any 
time t. The alternate or indirect role of impurities, e. g. in establish- 
ing initial densities of lattices vacancies, seemed unimportant when 
compared to the effect of stopping keV protons. The following premises 
thus seemed reasonable as a guide for initial experiments: 

(a) Proton damage restricted to a surface layer of the order of 
10“ 6 cm (100 A). 

(b) Measured radiative coloration sensitive only to selective 
impurities capable of functioning as recombination centers. 

(c) Measured radiative coloration dose- dependent, thus per- 
mitting laboratory simulation at a proton flux several orders 
of magnitude higher than solar. 










■$v 
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As a result of experimentation carried out during this second 
year the above premises require certain modification as follows: 

(a) Proton damage occurs in the bulk, probably as a result of 
e 1 - h* diffusion from a constant, p*- generated pair density 
in the surface. This bulk damage is the determining factor 
in radiative coloration and is slow to saturate. 

(b) The roles of impurities as (1) color centers and as 

(2) recombination centers affecting the kinetics of other 
intrinsic centers are not identical and must be defined, 

(c) Radiative coloration continues to appear dose-dependent only. 
The following discussion represents an attempt to illustrate and, where 
possible, conceptually formulate certain relevant optical damage 
kinetics induced by low energy proton irradiation. Evidence will be 
presented which illustrate charge carrier diffusion into the bulk followed 
by trapping at both native (intrinsic, e. g. , Schottky) and foreign 
(extrinsic) defects. 

4.-1 EVIDENCE FOR CHARGE CARRIER DIFFUSION 

In experiment (-71), a period of p -irradiation of a MgO single 
crystal resulted in decreasing the density of Fe from an initial value 
of 4. 1 x 10 ^/cm^ to a final value of ‘2. 3 x 10 ^/cm^ both values 

averaged over a measured bulk thickness »0. 1 cm. On the basis of 

* 2 2 
an incident surface of 1. cm , the bulk crystal volume 1 cm 7 a 

1 6 24 

volume containing 4. 1 x 10 Fe . These ions are initially uniformly 
distributed in the bulk. The volume associated with the p + - range 
( = 10 cm^) accounts for only 4. 1 x 10^ Fe"*" initially. Even if all of 
these Fe + in the "surface" layer changed valence state during the p + - 
bombardment, such changes, could not account for the measured 
decrease of 1. 8 x 10 16 fI + (= 4. 1 x 10 l6 - 2. 3 x 10 l6 ). 

Consistent evidence of this kind, involving net capture or release 

i • 

of an e or h by impurities located in the bulk as well as the surface, 
indicates e’ - h' pair diffusion into the bulk. Such diffusion results in 
the establishment of a significant density of free h* and e* at crystal 
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depths of the order of one-half or more of the crystal thickness x. In 
addition, data presented later indicates that such diffusion (as modified 
by e - h recombination and trapping) controls the growth rate of the 
observed color center bands. 

4. Z OPTICAL DAMAGE KINETICS 

A single damage model embracing all concurrent and sequential 

defect processes has proven to be unmanageable. This conclusion 

rests on the following basis: 

* » 

a. A solution of the e - h pair diffusion equation with recom- 
bination and trapping is very complex, 

b. Complete characterization of native defect structure is almost 
impossible. 

c. Many optical resonant states are not identified. * 

d. Precise effects of many substitutional impurities are not yet 
understood, e. g. AjK and Si in MgO. 

These complexities do not preclude limited analysis of certain observed 
electronic processes of real importance to the ultimate problem of 
achieving radiative stability in a white pigment. 

4, Z. 1 Band Gap Considerations 

An estimated energy level diagram is presented below (Figure 14) 
as an aid in following subsequent arguments. Note that only selected 
trapping and recombination levels are shown. These are established 
at thermal levels within the band gap based on available data or as 
inferred from reported qualitative thermal bleaching studies. Known 
thermal levels are: 

13 

V* = V, center : 1. 1 eV above valence band 
Mg 1 

. 4 

(OH - V Xir ) x = V^ TT center : 0.84 eV above valence band 
Mg OH 

x 14 

Vq s F center : 1. 53 eV below conduction band 

F + : <0.01 eV (i. e. unstable above 20°K)^ 

1 Q 

Note: See also Professor Kroger’s 10 electronic energy level schematic 
for MgO + Fe + Li based on Schottky disorder. 
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BAND 



k 3 =h -TRAPPING 
AND RECOMBINATION 

]f LEGEND: 

— Vq FILLED TRAP 

EMPTY TRAP 


BAND 


diagram for selected defects in MgO 



Hole trapping is designated by an arrow pointing upwards from 
the valence band, hole release by a reverse direction. Electron 
trapping is shown by an arrow pointing downwards from the conduction 
band, electron release by a reverse direction. In addition, the 
following symbols are used: 

H = 

H* = 

P = 

Ph = 

?H* = 
k = 

T = 

t = 

4. 2. 2 Growth Rate of Fe Band in the Bulk 

34 

Figure 15 illustrates the growth in the density of Fe in the bulk 

with p^ fluxes of the order of 3. x 10 ^ to 3. 4 x 10 ^ p + / cm^- sec. 

2+ 3+ 

This growth reflects a net conversion of Fe -*■ Fe whose rate at 
some particular point in the bulk may be expressed in terms of their 
densities as follows: 


total density of isolated with and without a trapped 

hole 

? i 

total density of OH - V ^ pairs, with and without a 
trapped hole 
density of free holes 

density of holes trapped at H = |Vj^J = centers! 

density of holes trapped at H""' = OH - 

[ v oh] centers 

2 — 

rate constant for trapping (cm /sec) = vc r 
time constant for thermal release from trap, sec 


time, sec 


d [Fe + ] _ , ft-2+1 

dt k J P(x) L Fe J (1) 

Estimation of the hole density depth- dependence p, as established by 

| * v“^ / 

e - h diffusion from a n surface layer 11 , constitutes the problem. 

This will be attempted parametrically, rather than rigorously, for 
reasons mentioned earlier. 
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IRRADIATION TIME t, 10 5 SEC 


3+ 

Figure 15. Measured growth of Fe density induced in the bulk of 
various MgO single crystals by 1. 0 keV proton irradiation 
(neutral beam) at a flux of 3. 0 x 10*- to 3. 4 x 10^ 

o'*'/ cm- sec 





(2) 


■ 

If U(x) represents a solution of the e - h pair diffusion 
equation, then varies as: 

P(x) = P S U (x) 

where: 

p g = hole density in surface layer 

Due to the high pair generation rate in the limited surface volume, at 
sometime, short compared to experimental times e. g. 6 x 10^ sec, 
the following occur: 

a* In the surface: 

(1) Defect densities saturate so that p - (p ) approaches 

suriace s 

a steady state value 
2+ 3 + 

(2) Fe /Fe equilibration occurs. Thus 




<x total Fe 


b. 


surface 

A fixed distribution of hole density from the surface layer 

into the bulk is established, i.e. U, , » constant. 

(x) 


Similar to earlier work we may write an expression for the 
time- dependence of p in the surface as follows: 


d P = 
dt 



k 2 p s 


(n v - p h *) + Pb*/t 


h 


- k.p 
1*8 


(H 


p h ) + p h /t i 


where: 

g 


i • 

proton generation rate of e - h in the surface 

19 3 

10 pairs/cm - sec 


(3) 


Note that upon early defect saturation in the surface layer (proton 
range) : dp/dt 0, P^j* — • * H‘, and p^ — - H. Based on reported 

thermal release values at room temperature of t, ~ 2. 5 x 10^ sec 

2 1 17 
and t 9 » 3. 5 x 10 sec and saturation values of p w * ~ 3. 5 x 10 / 

o ^ in 15 3 ^ 

cm' 5 and p„ * 0. 8 x 10 1 , p TT #/ t ~ 10 / cm - sec and p„/ t 

12T 3 ri c. ii 1 

« 3 x 10 1 /cm" 5 ~ sec. Note that these last two values are much 
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smaller than g. As a consequence, equation (3) approaches the 
following upon surface saturation: 


o “ g - P< 




p r « 


'~fI + ~~ 


Based on note a. (2) above, equation (5) may be restated as follows: 


(constant) g 
j^total F e 


n oc 
^S 


[total Fe] 


From equation (2) and the assumption = constant, we estimate that: 


(x) [total Fe] 


Rewriting equation (1) in terms of (7) yields: 


M. 


= k 3 C 


[ F * + ] 

total Fe 


and we observe that the growth rate of Fe in±he bulk (see Table 5) 

2 + 

depends on the ratio of the densities of Fe and total Fe. Since there 

is evidence (Figure 16) that in our experiments we are considerably 

removed from saturation in the bulk, we might compare experimental 

2 + 

rising kinetics on the basis of initial Fe densities in the bulk, i* e. 


48 



32 



0 20 40 60 80 100 
IRRADIATION* TIME t, 10 5 SEC 

3 + 

Figure 16. Estimated p + -induced growth to saturation of Fe in the 
bulk of two MgO single crystals 
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Table 5. Estimated -Fl* Growth Rate Dependence in the B.ulk 


Experiment 

No. 

[fI + L 

in 17 / 3 

10 /cm 

[tot. Fe] , 

17 3 

10 '/on 


■ 2+ 
Fe 

o 

[tot. Fe] 

(-66) 

42. 

42. 

1 . 


(-71) 

4. 1 

8. 7 

0. 47 

(-55) 

0. 7 

2. 7 

0. 26. 


Comparison of Figure 15 with the ratios in column four above 
indicates reasonable agreement between analysis and experiment a;s’ 
shown in Table 6, Thus we conclude that the rising kinetics for Fe 
approach linearity and are dependent on the ratio [fI + ] / [total Fe] . 


Table 6. Comparison Between Computed and Experimental 
Normalized Growth Rates for Fe in the Bulk 


Experiment No, 

d Jf1 + ] jdt 

(normalized to expt. -55 values*) 

Computed 

Experiment 

(-66) 

3.8 

2. 8 

(-71) 

1.8 

1. 7 

(-55) 

1 . 

1 . 


4. 2. 3 Growth Rate of V j Center 


In the bulk of the crystal, the following expres sion- for p^ ^ 
holds true: 


8p 


fo*) . = D - 8 -P _ k p 

3t n 2 R 3 p 

dx 


[fI + ] - k 2 p (H' 1 ' - p H *) + p H */ t 2 

“ kjP(H - p h ) + Pjj/Tj (9) 


where: p = p ( ^ 
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Previous assumptions (Section 4.2.2) of early surface saturation and 
early establishment of a constant p-distribution ranging from the sur- 
face layer into the bulk are retained. Thus as in equation (7), 

r v 1 

p (x ) 00 [total FeJ . Note that the diffusion term, first term in 

2 2 

equation (9), is positive since 3 p/9x (assumed constant) is always 
negative. The time- dependence of at a given x may be written as 
an ordinary differential as follows: 

d p H /dt = k !p( x ) (H - p H ) {10) 


where the back reaction, or thermal release, is neglected for rising 
kinetics. During the rising portion of the kinetics, except near the 
surface layer, it is likely that p^. « H and similarly Ppj* « H v . 
Therefore equation (10) may be approximated by: 


dp H /d‘ - kj p (x) H 


( 11 ) 


Substituting P( x j from equation (7) in equation (10) yields: 

1 


dp H /dt cc 


[total Fe] 


H 


( 12 ) 


It is important to note that H is crystal-dependent and is not 
influenced by irradiation. Such dependency is not directly measurable. 
Hence an alternate evaluation of H will be sought in V- radiation (in 
which case the dose is delivered to the crystal bulk ) carried to Vj and 
center satu ration . Time- dependent expressions for p, p^, and 
Ptt* thus uniformly distributed in the bulk at saturation are as follows: 


dp/ dt = 0 



- k 2 P [ H ‘ “ p H*] + p H* / t 2 

' k l p [ H ' p h] + p h /t i 


(13) 
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dp H /dt = kjp |h - p H - Pjj/Tj = 


= 0 


(14) 


dp H */dt = k 2 p [fT - p H «j - P H */ T 2 = 0 
Substitution of equations (14) and (15) into (13) yields: 


P =■ 


[ Fl+ ] a 


where: 


g = e - h pair generation rate 
subscript n indicates density at saturation 

Substitution of ( 16 ) into (14) followed by rearrangement yields: 


H = 


ML 

k, g 


k ,g 


vr 


n 


(P H ) 


n 


(15) 


(16) 


(17) 


Exact evaluation of equation (17) is prevented by a lack of values for 
and In lieu of this, it may be expressed as the proportionality: 


H ' [ Fi+ ] „ ( p n) n 


k l S 


N + ]n T 


(18) 


Combining equations (12) and (1-8) yields a single scaling expression 
for comparing experiments as follows: 


dp H / dt « 


1 


total Fej 


[ Fi+ J n (p H ) n 


k i g + _l 


Fe + 

L J n 


1 


(19) 
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On the basis of a Y-dose (absorbed) of approximately 1 Mrad/hr 
and an ionization requirement of ~20 eV per e ' - IT pair, pair genera- 
tion g is estimated as follows: 

10^ rad (100 erg/g)(6. 2 x 10^ eV/erg)(3. 58 g/cm^) 
hr (3600 sec/hr){20 eV/pair) 

15 3 

3. 1 x 10 pairs/cm -sec (20) 

\ 

Recall from earlier discussion that s 2,5 x 10^ sec. 

If one further assumes a value for k^/k^ ~ 250, the rising 

kinetics shown in Figure 17 maybe explained. Note that in the first 
12 

year’s effort , this ratio of k^/k^ was estimated as 1000. A value 
of 250 seems quite reasonable. 

Comparison of the center rising kinetics (dp^/dt) for experi- 
ments (-71) and (-66) are based on the following measured values: 

Experiment No. 

-66 
-71 

Table 7 which follows may then be constructed according to the scaling 
directed by equation (19). 

The result computed in column seven is that the product 
of scaling factors for experiments (-66) and (-71) are the same. This 
predicts that their rising kinetics for the center coincide, a pre- 
diction verified in Figure 17. The early induction period indicated is 
not yet understood. No saturation data for Y -irradiation of experiment 


10 Ions or Defects /cm 

[total Fe] _ F ® + __n_ V 1 q-- <PH>n 

42. 12. ‘ 1. 

8.7 1. (est. ) 0.79 


g “ 
g = 
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\D I AT I ON TIME t, 10 5 SEC 

3f Vl center induced in various MgO single 
irradiation (neutral beam) at a flux of 

3.4 by 1C)H p*/ cm^-sec 



Table 7. Scaling Factors for Comparing the Growth Rate of Vj 
Center (dV^/dt = dpfj/dt) in Experiments (-66) and (-71) 




Factors Which Reveal H 

dppj/ dt <x 

(Product 

Experiment 

1 


Curly 

Bracket 

in 

Equation 

(18) 


H oe 

No. 



N n 

Product 

of 

Columns 
3, 4, & 

5 

of 

Columns 
2 and 6 

- 66 

0. 238 

12. 

0. 06 1 

1. 0 

0. 731 

0. 18 

-71 . 

1. 15 

1 . 

0. 199 

0.79 

0. 157 

0. 18 


(-50) was available for correlation. It is interesting to note however 
that the early V} center generation rate observed in this case appears 
identical with those for (-66) and (-71), after a longer period of 'delay. 

4. 2. 4 Kinetics of Miscellaneous Centers 

The optical band peaking at 3.. 6 eV is thought to represent an 
intrinsic defect containing a trapped electron. Except for slightly 
different induction periods, early growth rates for this band appear 
crystal- independent (see Figure 18). Pending further identification 
of this center, no further correlation was attempted. 

4. 2. 5 Effect of Dose Rate on Generation Kinetics 

As discussed earlier, this work has involved exposure to proton 
dose rates of the order of 1000 X the solar wind proton flux. Time 
integrated fluences equivalent to an order of 30 years exposure in 
space were delivered. It is thus relevant to ask whether the optical 
damage would be equivalent between space and the accelerated 
laboratory experiments. 
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INCREASE IN MAX. ABSORPTION COEFFICIENT «, 


A EXPT. -71 (MODERATE PURITY-OXIDIZED) 
A EXPT. -66 (LOW PURITY-REDUCED) 

4 EXPT. -50 (LOW PURITY-AS RECEIVED) 


3.6 eV 

OPTICAL BAND 


IRRADIATION TIME t, l(r SEC 

Figure 18. Measured growth of 3. 6 eV band>in various MgO single crystals 
by 1. 0 keV proton irradiation (neutral beam) at a flux of 
3.0 x 10 11 to 3.4 x 10 11 p + /cm 2 -sec 





The argument will be related to the dependence of hole processes 
in the bulk, although a similar argument for electron processes could 
be made. The equation for hole processes in -the bulk is given again: 


9p (x, t) _ 

at 


-D 


d 2p (x, t) 
ax 2 



k 2 p(H* - p R *) - k lP (H- p H ) 


+ 




(9) 


In order to demonstrate strict dose dependence we must show that each 
term above is directly proportional to the hole generation rate g. Such 
demonstration permits g to be factored out as a simple time scale 
factor. 

In accord with earlier argument, after some short time period 
the distribution of free holes is given by: 


P 


(x) 


P u , > 

S (x) 


( 2 ) 


where: 



solution, of diffusion equation 


also 




( 21 ) 


For a given value of [total FeJ , note from equation (6) that: p g ce g. 
Therefore from equations (2) and (21) above we see that: 


9x< 


g 


(x) 
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Thus the first four terms in equation (9) above are shown to be 
proportional to g. When well below saturation, the thermal release 
terms (p^s/T^, P^/ T j) are of lesser importance and may be omitted. 
Since g is proportional to the proton dose rate g^ + , the factor g^ + is 
common to the four remaining terms in equation (9) and may be trans- 
posed as a time scale factor thus: 

8 p (*,.t> /S(t ‘V ) - 


We therefore conclude that these hole processes are independent of the 

dose rate g +. L.< 

p 

space simulation. 


dose rate Laboratory acceleration is thus valid in terms of 
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5.0 SUMMARY 


This work is associated with the development of improved 
thermal control coatings for space application. The effort and the 
conclusions derived therefrom maybe summarized as follows: 

a. MgO (both single crystals and polycrystalline material 
was selected as a typical dielectric solid transparent to 
solar (photon) radiation. Numerous specimens were 
exposed to a neutral., -1.0 keV proton beam in vacuum 
at a dose rate of 3.0x10^ to 3.4x10^ p + /cm^-sec 
(i.e. an order of 1000X the solar wind flux). Changes 
in optical and EPR spectra were measured in situ during 
exposure as a means of identifying proton -gene rated 
defects (important to radiative coloration) and determin- 
ing their growth kinetics. 

b. Material emphasis was placed on single crystals because 
of greater ease in the control of and identification of 
defect structure compared to powder s. Controlled physical 
processing (thermal) and chemical processing (e. g. oxida- 
tion, reduction, iron-doping) of specimens was-performed 
as an integral part of the experimentation. 

c. Resolution of proton -induced optical and EPR spectral 
displays has been carried out so as to identify important 
individual defect species and their densities. 

d. Measured changes in impurity defects establish the 
dominant coloration mechanism as consisting of e'-h pair 
diffusion into the crystal bulk from a surface layer (defined 
by the projected proton range) followed by recombination 
and trapping in the bulk. This conclusion is in contrast 

to the first year's theoretical effort wherein damage within 
the proton range was believed to determine radiative colora- 
tion. This diffusion-dependence greatly increases the 
influence of initial defect structure on the resultant optical 
damage. 
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e. Defect rising kinetics , i.e. well below saturation, are 
indicated to be dose -dependent. Thus laboratory 
simulation employing accelerated proton dose rates 
(relative to space) remain relevant to the work. 

f. Even on an idealized basis, a generalized model describing 
defect kinetics has not proven feasible. This results from: 

1. Lack of definition of many participating defect species. 

2. Difficulties in resolution and identification of defects. 

3. Complex solution to diffusion equation. 

g. Modeling by resort to scaling factors has proven both 
feasible and instructive. This has been used to correlate 

3-j. 

growth rate data for typical foreign defects (e.g. Fe ) 
and native defects (e. g. center) among various 
experiments. 

h. Comparison of proton-induced spectral absorption in 

MgO reveals a strong dependence on* initial purity (selective) 
and native defect structure. In this regard, the radiative 
stability of a processed, high purity crystal was 4 established 
by the small change, (*g) initial * 0.03 * 0.031, 

induced by a 2. 3 x 10 p + /cm^ fluence of 1 . 0 keV protons. 
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APPENDIX A 

ANALYSIS OF OPTICAL MEASUREMENTS 

1. 0 SINGLE CRYSTAL STUDIES 

Meaningful interpretation of proton-induced optical absorption or 
reflection changes (as indicated say in Figure 7} require resolution of 
the measured values into discrete component spectra, each indicative 
of a particular defect species- Further required is the establishment 
of a relation between measured, total spectral reflectance (R^h R) 
and changes in defect densities, e. g. by means of derived values of the 
total spectral absorption constant (n^~ or). 

1. 1 RELATION BETWEEN A* AND R 

Note that & sums all individual bands thus: 


“ = i [“ 


1 1- 4* etc 


■]■ [ 


r l N l + ^z^z + etc 


•] 


(A- 1 } 


The Fresnel reflection r at the air-surface interface is de_fin_ed in 
terms of the index of refraction n 1 of the crystal material: 


r = 


n -1 )/(n'-l-l}p. Optical ray schematics (Figure A-l) are useful for 
visualizing various terms in the relationship being sought. 

For wavelengths 0. 21<\<2. 5p, r<0. 10. As a consequence, 

2 

terms containing r are quite small and may be omitted. Terms in r 

and the large R term are as follows: 

& m 

Basis: (a) a is associated with total crystal thickness x. 

(b) Reflection at mirror second surface = R . First 
surface mirror reflection (Suprasil) is ignored. 

(c) Secondary back reflection resulting from multiple 

internal crystal reflection is accumulated as in 
Figure (b) above: Sr^ ” r 2 ^ r 2 + T^etc* 

(d) The attenuation factor e -crx represents the internal 
absorption in one pass normally through x. 


A-l 





r l = r 


(A2) 


I 


r 2 = 


r(l-r) 2 e- 2 “ x + [t( 1 -r) e’ 2 "*] r 2 ( 1 - r ) e' 2mt 


f , . . - 2a?xl 2 -Zax 

r{ 1-r) e r e 


2.. . -2 ax 

r ( 1 - r) e 


.2 -2Qfxf 1 2 -2ax , 4 -4ax 

r(l-r)e 1+re + r e 


-] 


(A3) 


Similarly: 


I 


^ -2ofx , t 2 -2orx t 4 -4cvx 

(1-r) e 1 + re +re 


m 


-] 


(A4) 


> , , ,4 2 - Zax f 

z—r . = r(l-r) R e 1 + 

4 m L 


2' - 2orx , 4 - 4o , x 

re + r e 


---•] (A5) 


I 


,, v6 ^2 -4cvx , . 2 -2ax , 4 -4ofx 
r_ - r(l“-r) R e 1 + re +re 

d m L 


--] 


(A-6) 


2 “ 2cvx 

The bracket factor is a geometric series with common ratio r e 

f 2 - 2^x1 " 1 

and first term = 1. The sum of such a series is: 1 - r e 

Summing (A2) through (A5) and substituting the series sum for the 
bracket factor yields: 


- L. 


r ~ r + 


1 M 

1-r e 


-2 ax r 

V -2»x -[ r(1 - r ) +R m (1 - r » 


_ , n . 6 „ 2 ”4orx 

4 2 I R e 

+ r(l-r) R 1 +* ~ = 

m I , 2 -2orx 


(A7) 


r e 
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For reasonable values of a, the third term in (A7) becomes small and 
can be dropped. Also the expression [l-r^ e -2cvx ] approaches a 
value of one. 


R ~ r + 


r(l-r) 2 + R m (l-r) 4 (l + rR m ) | e 


■] 


- 2 QfX 


(A8) 


Designating R q as the measured spectral reflectance before irradiation 
and R^ as that corresponding to some period of irradiation, equation 
(A8) may be written: 


R -r 
o 


[bracket factor] 


R^-r ~ [bracket factor] 
Dividing ( A9 ) by (A10) yields: 


-2a x 
o 


e 

-2a,x 

t 

e 


(A9) 

(A10) 


' (R 0 -r)/(R t -r) 


e 2 ( a t" a o'> K _ e 2A “ x 


This may be simplified to yield: 


Ah - — In 

2x 


R -r 
o 


R t~ r 


(All) 


(A12) 


From equation (A12) it can be seen that errors become 
significant at low. reflectance values (when R^r) and at high reflectance 
values (when R^R^). These represent the minority of cases. In the 
remainder errors of the order of 5 percent are considered typical. 
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Examination of measured reflectance values of unlr radiated high 
purity materials in. a spectral region where cr+o (e.g. \= 680 nm), 
indicated slightly low values (3 to 7 percent) as a result of losses in 
striations, between crystal joints in the mosaic array etc. Adding 
back these losses (on the assumption that they are wavelength- 
independent) and normalizing to with values calculated by (A7) 
when a - o) yields the estimated values of spectral absorption a ^ 
given in Figures 8 and 11. 

1. 2 GAUSSIAN ANALYSIS 

This analysis is based on the following set of spectral band 
parameters at 20°C : 


Band Width 


Band 

No. 

Designation 

Location of Band 
Peak, in eV 

at Half- Max. , 
W in eV 

Refere 

No. 

1 

Composite 

(unknown composition) 

5. 9 (est. ) 

0. 96 (est. ) 


2 

F + 

4.91 

0. 65 

7 

3 

c 2+ 

r 

4. 85 

0. 60 (est. ) 

15 

4 

f 3+ 

e 

4. 25 

0.99 

4 

5 

Unknown 

3. 6 

0. 60 (est. ) 

' 4 

6 

Unknown 

3. 1 

0. 44 (est. ) 

16 

7 

Unknown 

2. 7 

0. 52 (est. ) 

16 

8 

v i 

2. 3 

1. 07 

2 

9 

Unknown 

2. 2 (2. 15) 

0. 50 (est. ) 

4 


A major problem, as noted in the main text earlier, is the lack of 
definition in the near UV at (hv)>~5. eV. Several bands are believed 
to be in this region between 5 eV and the band edge of ~8 eV, including 
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3+ 

one known to be associated with an Fe transition. Additionally, the 

F and F bands are almost superimposed in the vicinity of 5.eV and 

2+ 

overlap the adjacent Cr band. As a consequence, spectral resolution 
in the region (hv)>4. 8 eV proved difficult and the results poor. 
Additionally, assumption of a Gaussian distribution may not be valid 
for certain of the observed bands. Estimated values indicated above 
were obtained by hand fit of experimental data. 

Let y = A a as obtained from (AZ2). Minimize the function, 

4> = ^ (y-y) 2 where y = measured value and y = the calculated value. 
Note that y can be expressed linearly in terms of each of the above 
nine components as follows: 


y = AjXj + A 2 X 2 --- etc. 


{ A13) 


whe re : 


X 



(E-E o ) 



(A14) 


c* - standard deviation 

' E = variable photon energy, eV 

E q = hv at component peak, eV 

Since the front factor in (A14) can be included in the coefficients A^, 
A^ etc. in (A13), we may write: 


X 


exp 


(E-E c r 

“ 17 " 


(A15) 


From the definition of a Gaussian distribution it can be shown that: 


— (A16) 

sfZ In 2 
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where: 


w = (E-E ) at half max. • 
v o' 

= W/2 

Computer solution of (A13) in terms of a least squares fit of the 
measured data was obtained. 
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APPENDIX B 


ANALYSIS OF EPR MEASUREMENTS* 


The physical basis for electron paramagnetic resonance depends 
on the spin magnetic moment or the orbital magnetic moment or both 
of an unpaired' electron about the nucleus of a paramagnetic species* 

In the presence of an external d-c magnetic field of strength H (in 
gauss), two quantum states will be associated with the combined mag- 
netic moment, p, of the unpaired electron. When electron spin moment 
is parallel to the field the energy is E = pH; when antiparallel, E = -pH. 
These two allowed states are thus separated by an energy difference 
AE = pH - (-pH) = 2pH. The value 2p may be replaced by the product 
gp, where g is a dimensionless variable called the spectroscopic 
splitting factor (and equals the ratio of the total magnetic moment to 
the spin magnetic moment) and p is a constant (=9. 27 x 10 ^ erg/ 
gauss) designated the Bohr magneton. The factor g is characteristic 
of the atomic environment of the electron and for a completely free 
electron g - Z. 0023. Larger values typically indicate trapped holes and 
smaller values trapped electrons. This splitting of electron energy 
levels by a magnetic field H may be shown schematically alFfollows: 



where 


AE = gpH 


(Bl) 


* Abstracted from: Blair, P. M. Jr. , Levin, H. , Honnold, 

V. R. , and Peffley, W. M. , "Study of the Combined Effects of Space 
Environmental Parameters on Space Vehicle Material, 11 George C. 
Marshall Space Flight Center NAS 8-21087 (Dec. 1969). 
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In order to effect a transition from lower to upper quantum energy 
level, electromagnetic radiation can only be absorbed at a frequency 
v given by: 


AE = hv 

_ 27 

where: h = Planck’s constant = 6* 63 x 10 erg-sec 
From the above one obtains the following: 

hv = g(3H 


(B2) 


and 


V = (gp/h)H 


(B3) 


For a free electron, the resonant or absorption frequency is: 


, = t2-°°Z3)9. . 27;x . lC5} , 2j 80H (B4) 

6. 63 X 10 

where: v is in MHz and H is in gauss 

Thus by use of a d-c magnetic field of 3400 gauss, one may obtain 
energy absorption at a convenient microwave frequency of the order 
of 9 500 MHz. 

The conventional procedure is to tune the microwave cavity con- 
taining the specimen to some precise frequency near 9500 MHz and then 
vary (scan) the magnetic field until absorption of the microwave energy 
is noted. The resonant signal thus obtained is displayed on a strip 
chart recorder, generally in the form of the derivative of the absorption 
line. This line can be studied in a way that is similar to methods used 
in other more familiar resonant phenomena. 
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Briefly, the following types of information on a radiation 
induced electron center can be deduced from a study of its resonance 
signal. If the resonance line is fairly narrow, the density of the 
centers can be estimated by a comparison of the intensity of the signal 
with that from a standard, usually an organic free radical. In a study 
of the kinetics of the formation and decay of a center formed by 
radiation, the change of the amplitude of the signal with time will give 
a direct indication of the dependence of the density with time. In some 
cases, the electron in a metastable state will interact with the nuclear 
spins of surrounding atoms. This interaction ( u hyperfine interaction 11 ) 
can cause a splitting of the observed resonance line. When this splitting 
is resolved, the number of hyperfine lines gives directly the nuclear 
spin of the nearby nuclei and, consequently, an indication of the position 
of the electron spin in the crystal structure. 

1.0 DETERMINATION OF g- FACTOR 

A knowledge of the g- factor, together with the width of the 
absorption line, leads directly to knowledge of the atomic structure 
surrounding the unpaired electron. Generally, a narrow resonance 
line (of the order of a few gauss) when associated with a g-factor 
close to that for the free electron is an indication that the electron is 
relatively isolated from the surrounding lattice. 

By determining the frequency at which the peak occurs, the 
g-factor can be determined to the degree of accuracy that v and H are 
known. From Equation (B3) an expression for g is 

g = (h/p)v/H 

= (6. 63 x 10“ 27 /9. 27 x 10" Zl )v/H 
= (0.7 1444) V /H 

where: v is in MHz and H is in gauss 
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2. 0 MEASUREMENT OF SPIN CONCENTRATION 

On the assumption that an unknown EPR line has been identified 
(at least tentatively) with a particular defect species, any understanding 
of its generation or annealing kinetics depends on a further knowledge 
of the defect concentration (density). Measurement of the number of 
spins in a given mass of a specimen of known density provides such 
data. 

The procedure observed in this laboratory was carried out in the 
order indicated below: 

a. A primary standard was selected with spin concentration, ' 
line width and location similar to resonance system of 
interest. 

b. A primary standard was used to calibrate a secondary stand- 
ard which was small enough to permit location within the 
cavity during EPR measurements on unknown specimens. 

c. With the calibrated secondary standard in place in the 
cavity, the unknown specimen is run in the cavity. 

Specifically a primary standard, consisting of powered pitch 

1 c 

dilated in powered KC1 to have 3 x 10 D spins, was obtained from 
Varian Associates. The secondary standard was a ruby crystal. It was 
cemented (styrene cement) in the microwave cavity at such an orienta- 
tion as to insure that its resonance line did not overlap the unknown 
lines displayed by the experimental specimen. 

The number of spins associated with the pitch in KC1 and ruby • 
lines may be compared on the basis of the area beneath their respective 
absorption lines . The EPR data is usually taken in the form of the first 
derivative of the absorption line. 
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Appropriate useof Equations (B6) and (B7) now permits integration of 
Equation {B5) by parts: 



A 


= uv - v du 


r°° co r® 

= / f(H-) dH =■ [f(H) . AH] ■- I AH - f'(H)dH 

J ^ — * CO 


(B8) 


Noting that’f(H) = 0 at H = co and H = -«>, Equation (B8) may be 
simplified: 


A = 



CO 

AH • f'(H) dH = 



AH • f T (H) dH + 



AH • f 1 (H) dHj 

(B9) 


Reference to the sketch, of f*(H) indicates that the product, AH * f'(H) 
is always negative* , since (either to the right or left of H Q ) when. AH is 
positive f 1 (H ) is negative and vice versa* Thus the integral in Equation 
(B9) is positive and the limits of integration are consistent with the 
designation of as the magnetic field intensity at the peak of the 
absorption line. 

Graphical integration under f^Hj.may be carried out by means of the 
app r oximatio n : 


2 AH - f'tHjAHi =2 (H.~ H o )f'(H)AH i (BIO) 

i i 

where the increment AHL is taken small and is not to be confused with 
the u lever arm 11 AH (=H> - H q ). 

Equations (B9) and (BIO) contain areas multiplied by n lever 
arms. n This fact points up the importance of obtaining narrow lines 
and/or good definition of line tails. Failure to do so ( as with broadened 
lines ) increases error rapidly, because of the magnifying effect of these 
lever arms in the integration. 
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The above procedure permits calculation of A(pitch) and 
A(ruby). The number of spins associated with the ruby crystal was 
then determined and used as a normalizing factor from which the 
number of spins associated with each of the unknown resonance lines 
in the experimental specimen was determined. 



